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ABSTRACT 

Extrapolation of perturbative nucleoli spin fluctuation rates seems to suggest a strong 
suppression of weak interactions in supernova cores. We derive a new sum rule for the 
dynamical spin-density structure function which relates the spin fluctuation rate to the 
average nuclear interaction energy. For a bremsstrahlung like structure function profile we 
show that instead of strongly decreasing, the neutrino scattering cross section is roughly 
density independent and axion emission rates increase somewhat slower than the lowest 
order emissivities towards the center of a hot supernova core. 



1 Introduction 



The cooling history of a newly born neutron star in the center of a supernova (SN) is 
mainly determined by neutrino diffusion. Numerical simulations employing the lowest order 
neutrino interaction rates calculated within the Glashow-Salam- Weinberg theory predict a 
cooling time scale which agrees remarkably well with the neutrino signal observed from SN 
1987 A [|]]]. The emission of novel weakly interacting particles like axions |2j could change the 
cooling time scale substantially which in turn allows to derive constraints on the properties 
of such particles ||. 

Within linear response theory weak interaction rates with a medium of nonrelativistic 
nucleons are determined, apart from the weak phase space, by only two dynamical structure 
functions, one for the density and one for the nucleon spin-density [|], [|. Some work has 
been devoted to their calculation but either the Landau theory of quasiparticles was applied 
assuming a "cold" nuclear medium 0, || or the authors focused on quasielastic scattering 
studying static structure functions (|, [7|, ||. Interactions of neutrinos and axions with a 
nonrelativistic nuclear medium are mainly governed by the local nucleon spin-density and 
its fluctuations. To lowest nontrivial order in the spin dependent nucleon-nucleon interactions 
causing these fluctuations, the relevant weak processes are of the nucleon bremsstrahlung 
type. Due to the Landau Pomeranchuk Migdal (LPM) effect @ which accounts for multiple 
nucleon scattering the inelasticity of these processes depends on the nucleon spin flip rate. 
In addition, once this rate becomes considerably larger than the medium temperature T, 
the total weak interaction rates tend to be suppressed [JTTJ. Since perturbative estimates 
for the nucleon spin flip rate can be as high as ~ 50T around nuclear densities, this could 
have profound implications for SN core physics || 10, 11, 12]. By dramatically reducing 
the predicted SN cooling time scale it would spoil the agreement between theory and the 
observed neutrino pulse from SN 1987 A JTTJ. On these phenomenological grounds it has 
been suggested that axial-vector neutrino scattering cross sections might be roughly density 
independent [0 instead of being suppressed at high densities by the LPM effect. 

In this letter we derive a new sum rule for the dynamical spin-density structure function 
(SSF) which provides an independent theoretical argument supporting this conjecture. It 
also predicts that emissivities for weakly interacting particles should increase somewhat 
slower than the lowest order rates at high densities. 



2 The Spin-Density Structure Function 

In terms of the nucleon field operator in the nonrelativistic limit, ip{x), the spin-density 
operator is given by cr{x) = \%1)\x)t , 4j{x) where r are the Pauli matrices. In the following 
we denote the momentum, coordinate, and spin operators for a single nucleon by p^, r^, 
and (7j, respectively, where i = 1, • • • , iV& runs over Nj, nucleons. Then, for a normalization 
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volume V, we can define the Fourier transform 



1 r 1 Nb 

a(t,k) = -j d 3 re- ik -V(t,r) = -^e'^a, 



In terms of these operators and the baryon density n& the SSF is defined as j5], 12 

4 r+oo 

S a (uj, k) = — / die 1 "' <<r(t, k) • <r(0, -k)) , (2) 

where (a?, k) is the four- momentum transfer to the medium. The expectation value (• • ■) in 
Eq. (0) is taken over a thermal ensemble. 

The contribution of S a to the neutrino scattering rate (per final state density) from four 
momentum (a>i,kx) to (a> 2 , k 2 ) can be written as \G\C\nb{?> — cos 9)S a (ui — co> 2 ,ki — k 2 ) 
with Gp the Fermi constant, Ca the relevant axial-vector charge, and 9 the angle between 
ki and k 2 . Similarly, the rate for pair production would read jG^C^n^ + cos 9)S a (— u)\ — 
u 2 , — ki — k 2 ) ||. The axion emission rate per volume, Q a , is governed by the same structure 
function [in an isotropic medium S^u, k) = S a (u, k) only depends on k = |k|]: 

® a = u N \2 b f2 I du^Svi-uj^) . (3) 

( 47r ) /a ^0 

Here, / a is the Peccei-Quinn scale and the numerical factor Cn depends on the specific axion 
model 0. Neutrino opacities and axion emissivities are therfore mainly determined by the 
SSF at thermal energies u ~ k <T '. 
Eq. (0) implies 



I 

J — i 



^ooSM k) = ~ ([H, <t(0, k)] ■ (7(0, -k)> , (4) 



where if is the Hamiltonian of the system of interacting nucleons for which we assume the 
following form: 

N b 2 i N b 

H = H + H int = Y / ^7 + ^T, V (*ij, <Ti, tr s ) . (5) 

1 = 1 iy^j 

Here, = r, — rj, M is the free nucleon mass, and V(rjj, <Tj, <Tj) is the spin dependent 
two nucleon interaction potential. For notational simplicity we restrict ourselves to only one 
nucleon species for the moment; the general case will be discussed further below. 

For free nucleons one gets f^(du>/27r) u>S a (u>, k) = k 2 /2M, in analogy to the well known 
f sum rule for the dynamical density structure function. In the latter case nucleon number 
conservation ensures that the f sum rule even holds in the presence of velocity independent 
interactions. In contrast, the f sum for the SSF is modified in the presence of spin dependent 
interactions since the nucleon spin is in general not conserved. 

For one nucleon species the most general two nucleon interaction potential is of the form 

El 

V(r, cr 1 , <r 2 ) =U(r) + U s (r)cr 1 ■ <r 2 + U T (r) (3cr 1 ■ f cr 2 ■ r - cr 1 ■ cr 2 ) , (6) 
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where r = 1*12, r = |r|, and r = r/r. We denote the spin dependent terms by V^j = 
Us{Tij)(Ti ■ o-j (the "scalar force") and V? = U^Tij) (3<Xj • <Xj ■ — cr^ • Cj) (the "tensor 
force"). In order to calculate the additional commutator in Eq. (^) from Eqs. ([I]), (13) and 
(|) we make use of the commutation relations erf, cr^ = i5ije abc a^, where z,j = 1 , - - • , iV& 
and e abc is the total antisymmetric tensor in the spatial indices a, b, c. After some algebra 
and using the symmetry properties of the Hamiltonian the modified sum rule reads 



£T %« s ^-m-m£ & + V S + « ■< ■ * - Kf) ) . ( 7 ) 



The kinetic nucleon recoil term is in general negligible compared to the ^/-dependent terms 
which govern the inelasticity of axial- vector interactions. If V(r, cr 1; <x 2 ) > —a/r s with a > 
and s < 2 the eigenvalues of H are bounded from below and the r.h.s. of Eq. ([?[) is finite 
as long as Us{r) and Ut{t) are integrable. This is the case for typical meson exchange 
potentials with hard core repulsion |TB|, |TJ], |T5|] . Assuming the three terms in Eq. (^) to be of 
similar size the r.h.s. of Eq. (^) is roughly proportional to the average interaction energy per 
nucleon W. At zero temperature and for SN core densities and compositions, W ~ 30 MeV 
corresponding to an average binding energy of about 10 MeV per nucleon. For T > 10 MeV 
nucleons are bound more weakly and W should be considerably smaller. We can therefore 
write 

— u>S a (co,li) ~AW< 100 MeV, (8) 

-00 Z7T 

where the inequality is a conservative bound reflecting our poor knowledge about the equa- 
tion of state for hot nuclear matter. Since it involves bound state energies, Eq. (§) is a 
nonperturbative result and will play an important role for the high density behavior of weak 
interaction rates below. 

The dependence on the momentum transfer k in Eq. (^) is expected to be only modest. 
In fact, for k < T < 50 MeV, we have |k ■ r| 1 within the range of the potential r s ~ l/m„ 
which is determined by the pion mass m n ~ 140 MeV. We can thus go to the long wavelength 
limit I, I, 0, m g, k -> 0, using S a (uj) = S a (uj, k -> 0). Eq. (0) then simplifies to 



+°° du 4 

— u;S a (u) = --(H T ), (9) 

where # T =§E$;^J- 

First, note that the scalar force does not contribute to Eq. ([|) because it conserves the 

total nucleon spin cr(0,k — > 0) [see Eq. ([[])]. Below nuclear densities the nucleon-nucleon 

(NN) interaction is dominated by one-pion exchange (OPE) leading to a tensor force. This 

contribution induces a spin orbit coupling and does therefore not conserve the total nucleon 

spin. Thus only the tensor force contributes to Eq. (HI) in the long wavelength limit. This 



agrees with the lowest order bremsstrahlung calculation for k — > [T5|]. Finally, note that 
the r.h.s. of Eq. @ is positive as it should be since the interaction induced correlations 
reduce (H T ) below the value for free nucleons, (H T ) = 0. 
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An additional sum rule ||, 
in the long wavelength limit: 

Note that for free nucleons Eqs. © and ([H]) yield S a (uj) = 2rr5(uj), whence only elastic 
scattering on the medium is possible in the absence of NN interactions. 



11, 12] can be obtained by integrating Eq. (H) and using Eq. 



3 Dilute Medium Limit 

At low densities, i.e. for large average inter-nucleon spacing, the interaction energy W in 
Eq. (H) is much smaller than the kinetic terms from the free Hamiltonian. In case of the 
long wavelength limit, Eq. ([]), we can therefore treat Ht as a small perturbation and write 
to lowest non-trivial order in Ht- 



(11) 



Here, E®, \n) are the eigenvalues and eigenstates of the free Hamiltonian Hq, respectively, 
\n) 1 are the eigenvectors of H + Ht to first order in Ht, and Z = J2n ex P ( — E^/T) is 
the normalization factor. Assuming nondegenerate eigenstates for simplicity and applying 
standard first order perturbation theory for \n) 1 we can express everything in terms of zeroth 
order quantities. Dropping the index from now on, Eq. ([TT]) reduces to the negative definite 
expression 

I P -En/T _ —E m /T 

(Ht) = -= E £ — p 1 K^tUI 2 , (12) 

where (H^mn =o (^| Ht \n) . This matrix element is expected to vary in E m — E n over a 
scale > 3m^/M ~ 50MeV where m n is a typical momentum scale in the NN interaction 
potential. Therefore, for T < 50 MeV the thermal factor in Eq. ( |12|) can be approximated 
by 5(E m — E n ). Converting the sum over m into an integral over E m Fermi's golden rule 
finally gives W ~ —(H T )/N b = Y a /{2^). Here, r CT is the average perturbative NN scattering 
rate mediated by Ht which is a measure for the spin fluctuation rate. The spins fluctuate 
on a time scale given by the inverse energy scale of the tensor force which causes the spin 
fluctuations. At high densities we use IV = — 2ti{Ht) /Nf, as an effective spin flip rate. 



4 Saturation of Spin Fluctuation Rates 

We now use the sum rules Eqs. (^) and ( |T0"D to determine the qualitative form of S a (u>, k) in 
the long wavelength limit. To this end let us introduce the dimensionless quantity S a (x) = 
TS a (xT) with x = lj/T as in Ref. fill. Due to the principle of detailed balance, S a (uj, k) = 
S a (— to, — k)e W//T , it is sufficient to specify S a (x) for x > only. 
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Introducing the dimensionless effective spin flip rate 7 CT = T a /T, we can write the sum 
rule Eq. @ as 

f+°° cb ~ / 2 7(T AW 

/ — (1-e J = -^r, (13) 

JO v ' 7T -/ 

where in a newly born neutron star 7^ does not increase beyond a few. 

Furthermore, since in a hot SN core the thermal energies are expected to be considerably 
higher than the interaction energy W, within a first approximation we can neglect spin 
correlations in the second sum rule Eq. (|10J) and write 

r + OO Ann . 

/ ps a (x) (l + e-) ~ 1 . (14) 
Jo zi v ' 

For the following discussion we consider the general case of an ensemble of neutrons 
and protons with fractional number densities Y n and Y p . Introducing the isospin operators 
Ti for nucleon i, cri in the definition of S a [see Eqs. @ and (0)] has to be multiplied by 
[1 + (tj) 3 ] C j 4 iP /2+[1 — (Tj) 3 ] Cj4,n/2. Here, Ca, p and C^n are the relevant proton and neutron 
axial-vector charges. Moreover, there will be additional terms proportional to Tj • Tj in the 
interaction potential Eq. (^). However, this leaves our discussion qualitatively unchanged 
since the additional isospin operators appearing under the expectation values only lead to 
additional factors of order unity. If correlations among different nucleons are absent the 
r.h.s. of the sum rules Eqs. (0) and ([14]) get multiplied by (Y p C\ p + Y n C\ n ) / (C 2 A;p + C\ n ). 

Parametrizing the high u behavior of S„ by S a {uj) oc u^ n , classical collisions would lead 
to n — 2. On the quantum mechanical level the deviation of S a (u) from 2tt8(u) is to lowest 
order in the strong interactions given by nucleon bremsstrahlung. Using a dipole like OPE 
potential without a hard core cutoff yields n = 5/2 and n = 3/2 in the case of one and two 
nucleon species, respectively |5|, [l|, 0, |i~7| . The non-existence of the f sum Eq. flilf ) in the 
latter case stems from the unphysical r -3 divergence of this potential at r = 0. Except for s 
waves this divergence indeed leads to an infinite (Ht)- If one regularizes the potential by a 
hard core repulsion f sum integrability is restored. 

This motivates the following representative ansatz: 

S a {x) = -jj^ forx>0, (15) 



where a and b are positive constants. The sum rule Eq. (|T3| ) is sensitive to the high energy 
behavior and therefore mainly to a. In contrast, the sum rule Eq. ([14]) probes the "infrared" 
regime which is sensitive to h. Eq. (|T^) is of the form expected from nucleon bremsstrahlung 
where b accounts for the LPM effect. 

We can now pick a number for a, determine the corresponding value of b numerically from 
Eq. ( |i~4l) and compute the f sum Eq. ( fL3|) . The result is plotted in Fig. 1 as a function of a. 
Most importantly, from the expected density and temperature dependence of W we expect 
the f sum to increase monotonically towards the SN core before saturating at a value of 
order unity. As a consequence, the thermally averaged axial- vector neutrino scattering cross 
section {a a) which dominates the neutrino opacity should roughly scale as T 2 being density 
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independent as naively expected (see Fig. 1). Furthermore, the axion emission rate from 
Eq. (HI) approximately scales as nf } T a T i . The lowest order axion emissivities should therefore 
be multiplied by Tv/T^ whenever this ratio is smaller than 1. Here, is the lowest order 

1 /2 

spin flip rate extrapolated from the dilute medium limit. For example, T a ~ 32 MeVpi4T 1( 



10 

114 , 



for the standard OPE calculations |fL5[ |16|, [17| , where pu is the mass density in 10 gem" 
and T 10 = T/lOMeV. A turn over in (<ja)/T 2 and Q a /(nbT 3 ) typically only occurs at 
7 CT > 10 and is the less pronounced the stronger S a (uj) falls off at large u. The absence of a 
decrease of these quantities at high density is therefore rather independent of uncertainties 
in the exact saturation value for r y CT . 



5 Summary 

Neutrino opacities and axion emissivities are governed mainly by the SSF. We have derived a 
new sum rule for the SSF which corresponds to the f sum rule for the density structure func- 
tion but depends on the nucleon spin flip interactions. Our treatment so far assumes absence 
of possible pion and kaon condensates. Employing an infrared regularized bremsstrahlung 
spectrum for the functional form of the SSF we have shown that the effective spin fluctuation 
rate T a must saturate somewhere below ~ 150 MeV which is within factors of a few of SN 
core temperatures. Neutrino scattering cross sections should therefore exhibit the naive T 2 
scaling whereas axion emissivities should increase somewhat slower than the lowest order 
rates at high densities. There is no turnover of weak interaction rates towards the SN core. 
These results have an important impact on SN cooling simulations and their application to 
the derivation of axion mass bounds. They are also relevant for the rates for URCA processes 
and emission of right handed neutrinos. 
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Figure Captions 



Fig. 1. The f sum Eq. (|T3| ) characterized by 7 CT as a function of the parameter a (solid line). 
Also shown in arbitrary units are the axion emission rate per baryon Q a /nb (dashed line) 
and the thermal axial-vector neutrino scattering cross section {a a) (dotted line) normalized 
to a fixed temperature. The physical range is where the f sum is smaller than a few. 
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